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We present a detailed formalism of the microscopic particle-rotor model for hypernuclear low-lying 
states based on a covariant density functional theory. In this method, the hypernuclear states are 
constructed by coupling a hyperon to low-lying states of the core nucleus, which are described by the 
generator coordinate method (GCM) with the particle number and angular momentum projections. 

We apply this method to study in detail the low-lying spectrum of iC and ^aNc hypernuclei. We 
also briefly discuss the structure of ^^ISm as an example of heavy deformed hypernuclei. It is 
shown that the low-lying excitation spectrum with positive parity states of the hypernuclei, which 
are dominated by A hyperon in s-orbital coupled to the core states, are similar to that for the 
corresponding core states, while the electric quadrupole transition strength, B{E2), from the 2f 
state to the ground state is reduced according to the mass number of the hypernuclei. Our study 
indicates that the energy splitting between the first 1/2“ and 3/2“ hypernuclear states is generally 
small for all the hypernuclei which we study. However, their configurations depend much on the 
properties of a core nucleus, in particular on the sign of deformation parameter. That is, the first 
1/2“ and 3/2“ states in are dominated by a single configuration with A particle in the p-wave 
orbits and thus providing good candidates for a study of the A spin-orbit splitting. On the other 
hand, those states in the other hypernuclei exhibit a large configuration mixing and thus their energy 
difference cannot be interpreted as the spin-orbit splitting for the p-orbits. 

PACS numbers: 21.80.-|-a, 23.20.-g, 21.60.Jz,21.10.-k 


I. INTRODUCTION 

The development in A-hypernuclear spectroscopy has 
enabled one to explore several aspects of hypernuclear 
structure [1, 2]. Since hyperon-nucleon and hyperon- 
hyperon scattering experiments are difficult to perform, 
the information on the AN and AA interactions have been 
extracted from such studies. Moreover, since a A hyperon 
is free from the Pauli principle from other nucleons, it 
can go deeply inside a nucleus, which can also be used 
as a sensitive probe in order to study the structure of 
normal nuclei. Theoretically, many methods have been 
developed to investigate the spectroscopy of hypernuclei, 
such as the cluster model [3-9], the shell model [10-12], 
the ab-initio method [13], the antisymmetrized molecular 
dynamics (AMD) [14-17], and self-consistent mean-field 
models [18-26]. Among them, the self-consistent mean- 
field approach is the only microscopic method which can 
be globally applied from light to heavy hypernuclei. 

One of the characteristic features of atomic nuclei is 
deformation in the body-fixed frame. In the mean-field 
approach, the optimized deformation is automatically ob¬ 
tained by minimizing the total energy of a system in the 
mean-field approximation. It was shown, however, that 
the potential energy surface of a hypernucleus is generally 
softer against deformation than that of the corresponding 
core nucleus [21]. This implies that the shape fluctuation 
effect, which is not included in the pure mean-field ap¬ 


proximation, will be more important in hypernuclei than 
in normal nuclei. Furthermore, in order to connect mean- 
field results to spectroscopic observables, such as B{E2) 
values, one has to transform the results to the laboratory 
frame. To this end, one has to rely on additional assump¬ 
tions such as the rigid rotor model, which however would 
not work for, e.g., nuclei with small deformation or with 
shape coexistence. To quantify the impurity effect of A 
particle on nuclear structure, one thus has to go beyond 
the pure mean-field approximation. 

In our previous publication [27], we have proposed a 
new approach using a microscopic particle-rotor model 
(PRM) for the low-lying states of single-A hypernuclei. In 
this method, the A particle is coupled to the core nucleus 
states while the A hyperon interacts with the nucleons 
inside the nuclear core. For the core nucleus, a beyond- 
relativistic-mean-field approach is applied for low-lying 
states by carrying out the angular momentum and the 
particle number projections as well as the configuration 
mixing with the generator coordinate method (GCM). 
We have successfully applied this method to the spectrum 
of A^e. 

The motivation of the present work is to introduce a 
detailed formalism of this method. At the same time, we 
also apply it systematically in order to study the low- 
lying states of single-A hypernuclei in the mass region 
from light to heavy. For this purpose, we will first dis¬ 
cuss the hypernucleus, which is an ideal hypernu¬ 
cleus in order to discuss the spin-orbit splitting. We will 
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FIG. 1: (Color online) A schematic picture of the microscopic 
particle-rotor model for A hypernucleus, in which r denotes 
the coordinate of the A hyperon. In this approach, the nu¬ 
clear core states are described microscopically with the multi¬ 
reference density functional theory. 


then discuss example of hypernuclei in the 

sd-shell region with a prolate deformation. The ^®^Sm 
hypernucleus, which is a well-deformed system and has a 
well-developed ground-state rotational band, will also be 
considered as an example of heavy hypernuclei. 

The paper is organized as follows. In Sec. II we present 
a detailed introduction to the formalism of microscopic 
PRM based on the covariant density functional theory 
(CDFT). In Sec. Ill, the results for the low-lying states 
of hypernuclei and the corresponding core nuclei will be 
presented. Finally, we will summarize the conclusions of 
this paper in Sec. IV. 


II. MICROSCOPIC PARTICLE-ROTOR MODEL 
FOR A HYPERNUCLEI 


with r and being the coordinates of the A hyperon 
and the nucleons, respectively. J is the angular mo¬ 
mentum for the whole system while M is its projection 
onto the z-axis. is the spin-angular wave func¬ 

tion for the A hyperon. |$/„) is the wave functions of 
the low-lying states of nuclear core, where I represents 
the angular momentum of the core state and n = 1 , 2 ,... 
distinguish different core states with the same angular 
momentum /. For convenience, hereafter we introduce 
the shorthand notation k = {jUn} to represent different 
channels. 

In the relativistic approach, ^k(j') is the radial wave 
function of a four-component Dirac spinor and it can be 
written in the following form 


^k{r) 


( fk{r) '\ 

\igk{r)cr -rj ' 


( 3 ) 


The Hamiltonian H for the whole A hypernucleus can 
be written as 


(NA), 


H = fA + f^ 

2=1 


Hr 


( 4 ) 


where Ac is the mass number of the core nucleus. The 
first term in Eq. (4), Ta, is the kinetic energy of A hy¬ 
peron. 


Ta =-ia • Va -I-7 °mA. (5) 

Here, m\ is the mass of the A hyperon and a and 7 ° are 
Dirac matrices. The second term in Eq. (4) represents 
the NA interaction term between the valence A and the 
nucleons in the core nucleus. It is composed of both a re¬ 
pulsive vector-type term and an attractive scalar- 

type term , for which we take the following contact 

coupling forms. 


A. Coupled-channels equations 

The basic idea of the microscopic PRM for a single- 
A hypernucleus is that the valence A hyperon couples 
to the low-lying states of nuclear core, as is illustrated 
in a schematic picture of Fig. 1. In this approach, a 
hypernucleus is described in the laboratory frame and the 
wave function of the whole A hypernucleus is constructed 
as 

'i'jM{r,{rN}) = ( 1 ) 

where 

Wl) = (2) 


V^^^\r,rN) = ay^6{r - tn) ( 6 ) 

V^^^\r,rN) = - rNh%. (7) 

These NA interactions correspond to the leading-order 
four-fermion coupling terms in the effective interaction 
proposed in Ref. [28] . For simplicity, the possible higher- 
order derivative and tensor coupling terms are neglected 
in the present study. We note here that the spin-orbit 
interaction of hyperon is automatically taken into ac¬ 
count in the relativistic framework without introducing 
an additional parameter. The last term in Eq. (4) is 
the many-body Hamiltonian for the core nucleus satisfy¬ 
ing the equation i7c|4>/„) = In this paper, we 

use He that corresponds to the following point-coupling 
energy density functional (EDF) [29], 
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£'c[{pj,{jf}] = Tr[(Q; • p +/3m)py] + J 


4 + |p| + f 4 + 'fpsAps + f 

+ ^^^jTvijTv)tJ. + ^^jTV^UTv)f, + ) ( 8 ) 


where is the four-component electromagnetic field, 
and the densities pi and currents are bi-linear combi¬ 
nations of Dirac spinors, namely with i = S,V, TV 
representing the symmetry of the coupling. The sub¬ 
script S stands for isoscalar-scalar (Tg = 1), V for 
isoscalar-vector (Fy = and TV for isovector-vector 
(r^y = type of coupling characterized by their 


transformation properties in isospin and in space-time. 

In order to obtain the radial wave function given by 
Eq. (3) and the energy of hypernuclear low-lying states, 
we multiply to the total Schrddinger equation, 

= TjI’I'jm}, and integrate it over f and {rjv}. 
This leads to the following coupled-channels equations. 


d K — 1 
dr r 


g,{r) + - Ej)Mr) + ^ [c/f'(r) + Uf{r)] f,,{r) = 0 


- Ej)gk{r) - {r) gk'{r) = 0 


(9a) 


(9b) 


where n is defined as k = + 1/2). With 

the multipole expansion for the S function in coordinate 
space, 

S{r - r'} = ^ (10) 

X,/j. 

the vector and scalar coupling potentials in Eqs.(9a) and 
(9b) read 


uY(r) ^ 




J I j 
A f I' 


NA i-ni-r,' P \ 

x«v P\V i'^)^ 


( 11 ) 


and 


Uf{r) ^ 


2=1 


= (-i)^'+'+"Eo'^iii^^ii/^') 


JI J 
A f I' 


NA Inl„> / \ 

x^s Pa.s 


( 12 ) 


scalar transition densities defined, respectively, as 


p 1 E '(0 = (' J’/JiE 


2=1 

Ac 


S{r - TMi) 

TNiT 


U(r7v0||$7 ,),(13a) 


p1 E'(0 = 11 E EA(f iv.)| |$7„a3b) 


i=l 


TNiT 


between the nuclear initial state Ini and the final state 
In ■ The detailed expression for the transition densities in 
the non-relativistic multi-reference DFT framework has 
been derived in Ref. [30] by one of the present authors. 
The formalism has also been generalized to the relativis¬ 
tic case within a multi-reference CDFT (MR-CDFT) to 
study a “bubble” structure in light nuclei [31, 32]. In this 
work, we extend this formalism to study low-lying states 
of hypernuclei. 

With the radial wave function I%k{r) in the coupled- 
channels equations (9a) and (9b), one can compute the 
probability Pk of the channel k in the total wave function 
^'jM as, 

Pk = J r'^dr\^k{r)y 
= J r^dr [\fk{r)y + Igkir)^] ■ 


where pE"' P^\Y (^) reduced vector and The wave function is normalized as Pk = 1. 


Inl 


( 14 ) 
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B. Projected potential energy surface 

In order to apply the formalism presented in the pre¬ 
vious subsection, one has to specify the core states . 
A simple choice for this is to construct them as the pro¬ 
jected mean-held states with the intrinsic deformation /3 

\^IMrm=PkKP^P^\m), ( 15 ) 

where the particle number projector has the form, 

1 /■27r 

P^-=— / , {T=n,p) (16) 

27r Jo 

and the operator P^^x i® three-dimensional angular 
momentum projection operator given by 

PmiK = g^2 J (17) 

Here, H represents a set of Euler angles {(j),9,ijj), and 
the measure is dH = dcj) sin 9d9dtp. P(H) and 
are the rotation operator and the Wigner D-function, 
respectively [33]. In Eq. (15), the wave function 
|(^(/3)) is a Slater determinant of quasi-particle states 
with quadrupole deformation /3 generated with the con¬ 
strained relativistic mean-held (RMF) calculation. For 
simplicity, in this paper we consider only the axial de¬ 
formation for the nuclear core and thus the K quantum 
number is zero in Eq. (15). 

In this approach, the hypernuclear states '^jm are 
given for each deformation fj. This allows one to con¬ 
struct the projected energy surface for hypernuclei. That 
is, one can dehne the total energy Ej(/?) by taking 
{'i' jm\H\'^jm) at each deformation /3 and for each spin- 
parity, . Notice that the coupled-channels equations 
(9a) and (9b) are solved at each deformation, and thus 
the effect of core excitations is taken into account in the 
projected energy surface so obtained. 


C. Multi-reference covariant density functional for 
nuclear core states and transition densities 

One can improve the calculations in the previous sub¬ 
section by using from the MR-CDFT calculation in 
the context of generator coordinate method (GCM) [34- 
36], that is, 

Wi.M,) = E F^.mPL.KP'^p^wm) ( 18 ) 

/3 


to an arbitrary variation of which leads to the 

Hill-Wheeler-Griffin (HWG) equation [37], 

^ [^"(/3,/3') - Pz„^"(/3,/3')] F^iP') = 0. (19) 

P' 

Here, ^"(/3,/3') {^mP^oP''P^ and 

= {(p{P)\PlPoQp^P^\(p{P')) are the norm 
and the energy kernels, respectively. In the calculations, 
the energy overlap in the energy kernel is taken to be 
the same functional form as in the nuclear mean-field 
energy, but replacing the densities and currents with 
mixed ones, that is, off-diagonal components of the 
density and current matrices [34-36]. 

Since the projected mean-field states do not form an 
orthogonal basis and the weights Fj^{P) in Eq. (18) are 
not orthogonal functions, it is convenient to construct a 
set of orthonormal collective wave functions giP as [38] 

aim = E (20) 

p' 

Notice that the modulus square of gi{P) does not rep¬ 
resent the probability to find the deformation /3 in the 
GCM state. For the axial symmetric case, however, gi{P) 
provides a good indication about the dominant configu¬ 
rations in the collective states. 

With the GCM wave functions for one can cal¬ 

culate the reduced vector transition density in Eq. (13a) 
as follows [30], 

^ P,P’ 

X E(70Aj^|/'0) 

Ky 

X /dfr;,(f)((^(/3)|py(r)Po4P^P^|<p(/3')), 

( 21 ) 

where the notation I = y/2I -|- 1 is introduced for simplic¬ 
ity, and the vector density operator is defined as follows 

Ac 

py(r) = E'^(i’- ( 22 a) 

i=l 

The scalar reduced density in Eq. (13b) can also be ex¬ 
pressed in a similar way. 

D. Electric quadrupole transition strengths 
between hypernuclear states 


In this wave function, a set of Slater determinants with 
different quadrupole deformation P is superposed accord¬ 
ing to the idea of GCM. The weight function FPI^[P) in 
the wave function (18) is determined by requiring that 
the energy expectation value is stationary with respect 


The electric quadrupole {E2) transition strength from 
an initial state \ Ji) to a hnal state | J/) in A hypernuclei 
is defined as 


B{E2]J, Jf) 


1 


2A-kl 


(J/||Q2||<A) 


(23) 
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Here, the E2 operator reads Q 2 ti = Sub- 

iGp 

stituting the wave function for the hypernuclear states, 
Eq. (1), to this equation, one finds the reduced matrix 
element to be 

ki,kf 

(24) 

with 

( 25 ; 

(see Eq. (7.1.8) in Ref. [33]). Here, (/ny:||Q 2 ||-fni) is 
the reduced E2 transition matrix element between the 
nuclear core states |//,n/) and and it is related 

to the proton vector transition density in Eq.(21) as 

{Inf\\Q2\\Ini) = ii J drr‘^pl"^'{r). (26) 

III. RESULTS AND DISCUSSIONS 

We now numerically solve the coupled-channels equa¬ 
tions and discuss low-lying spectrum of hypernuclei. The 
procedure of the calculations and the numerical details 
are listed as follows. 

(i) Self-consistent deformation constrained RMF-hBCS 

calculation for the nuclear core states: This step is to 
generate a set of deformed states \^p{l3)) with different 
quadrupole deformation /?. The Dirac equation for nu¬ 
cleons is solved with the basis of a three-dimensional har¬ 
monic oscillator (3DHO) with = 10 major shells. 
The oscillator length parameter in the 3DHO is cho¬ 
sen as bx = by = bz = y/hJrrujjQ, where m is the nu¬ 
cleon mass and the oscillator frequency is determined 
to be hujQ = MeV. In the most of calculations 

shown below, we employ the non-linear point-coupling 
EDF with the PC-Fl [29] set for the particle-hole chan¬ 
nel, although we also use the PC-PKl set [39] for ^'’Ne 
in order to study the parameter set dependence. In these 
energy density functionals, a density independent S force 
is used for the particle-particle channel, supplemented 
with an energy-dependent cutoff for the pairing active 
space [40]. 

(ii) MR-CDFT calculation for the low-lying states of 
nuclear core: This step is to obtain the wave functions 
<!)/„ for the core state In- The mean-field wave functions 
are projected onto good particle numbers (N, Z) and an¬ 
gular momentum /, which form a set of non-orthogonal 
basis. The Gauss-Legendre quadrature is used for inte¬ 
grals over the Euler angle 0 in the calculations of the 
norm and Hamiltonian kernels. For the nucleus, 


the number of mesh points in the interval [ 0 , tt] for 9 
is chosen to be 14. The number of gauge angle p for 
the particle number projection is chosen to be 7. For 
^°Ne and we use the mesh points of 16 and 9 

for the anglar momentum and the particle number pro¬ 
jections, respectively. The energy and wave function of 
the core state are determined by solving the HWG 
equation, Eq. (19) [34-38]. With the wave functions of 
nuclear core states, one can calculate the transition den¬ 
sities which are used to determine the coupling potentials 
in the coupled-channels equations. 

(Hi) Coupled-channels calculation for the low-lying 
states of A hypernuclei: With the coupling potentials 
so obtained, the coupled-channels equations are solved 
by expanding the radial wave function (r) on the 

basis of eigenfunctions of a spherical harmonic oscillator 
with 18 major shells. From the solutions of the coupled- 
channels equations, we construct the spectrum of hyper¬ 
nucleus and calculate the B{E2) transition strengths. 

A. Application to 

1. Properties of the core states 

We first apply the method to ^^G. Figure 2 shows the 
energy curves for the core nucleus, The dotted line 
is obtained in the mean-field approximation, while the 
other lines show the projected energy surface for each 
angular momentum /. The mean-field energy curve ex¬ 
hibits a pronounced minimum at the spherical configu¬ 
ration with a steep rise with deformation /3, as expected 
for a nucleus with large neutron and proton shell gaps. 
The energy gained from restoration of rotational symme¬ 
try increases with deformation (3 and together with the 
particle number projection the location of energy mini¬ 
mum is shifted on the curve. The minimum of the energy 
curve with O'*" is found on the oblate side with (3 = —0.3. 
Besides, the second minimum appears around (3 = 2.4 in 
the mean-field energy curve, which is shifted to (3 = 2.7 
in the projected energy curve for /’^ = 0+. It has been 
shown that the configuration for this minimum has a 3a 
linear-chain structure [41]. 

After mixing the projected mean-field configurations, 
one obtains the energies of low-lying states (see the filled 
squares in Fig. 2). The wave function of these states 
is displayed in Fig. 3 for the lowest three states with 
/ = 0,2,4, and 6 . The ground state of is domi¬ 
nated by the spherical configuration. The collective wave 
functions and the energy spectrum indicate that there is 
a coexistence of an anharmonic spherical vibrator and 
an oblate deformed band at low excitation energies of 
^^G. Both structures are not pure and distorted by their 
strong mixing. The high-lying O 3 ", 2^ and 4^ states seem 
to form a rotational band dominated by the 3a-linear 
configuration, in which the collective wave functions are 
much extended to a large deformation region. Similar ro¬ 
tational band corresponding to a 4a-linear configuration 
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FIG. 2: (Color online) The energy curve of the mean-field 
state and the particle-number and angular momentum pro¬ 
jected = O'*', 2"*', and 4^) states for as a function of 
the intrinsic quadrupole deformation (3. The filled squares 
indicate the lowest three GCM solutions for each 7'^, which 
are plotted at their average deformation ^ |Sn(/3)P/3i 

where gh{l3) is the collective wave functions defined by Eq. 
(20). The results are obtained with the PC-Fl force. 



FIG. 3: (Color online) The collective wave functions gi, given 
by Eq.(20), for the first three states in with spin-parity 
of 0+,2+,4+ and 6+. 


has also been found in the high-lying states of [42]. 

Figure 4 shows the vector and scalar transition densi¬ 
ties in the low-lying yrast states (n = 1) of where 
the multipolarity A is taken as 0,2, and 4. Po°(r) is noth¬ 
ing but the total nucleon density for the 0]*' ground state 
multiplied by a factor It is shown that the transi¬ 

tion density decreases by one order-of-magnitude as 
A increases from 0 to 2, and from 2 to 4. Besides, we 
also plot the transition densities with A = 0, 2, and 
4, with A = 2,4, and 6, and with A = 4, 6, and 8. 
Notice that the vector and scalar transition densities are 


slightly different from one another. 

The convolution of the proton vector transition densi¬ 
ties in Fig. 4 with a Gaussian form factor for a finite pro¬ 
ton size yields the charge transition densities, P^j^^h 
which are related to the form factor F^lq) for electron 
scattering with an angular momentum transfer L by the 
following relation [30], 

/a ' 

= drr^PL:k'ir)jL{qr), (27) 

where jhiQT) is the spherical Bessel function. The co¬ 
efficient V^/Z is chosen so that the elastic part of the 
form factor Fq^q) is unity at g = 0, as in the experimen¬ 
tal data. A comparison of our results with the data is 
shown in Fig. 5. One can see that the form factors 
Fhiq) are in rather good agreement with the data ex¬ 
cept for the underestimation of the elastic form factor 
after the first minimum, as was found also in the recent 
studies for [43] and ^^Mg [30] based on the Skyrme 
forces. This may be because the spreading of the col¬ 
lective wave function in quadrupole deformation space is 
somewhat overestimated in the calculations, decreasing 
the weights of the large-g components of the transition 
density [30, 43]. In fact, the charge radius of by the 
present GCM calculation is 2.57 fm, which is larger than 
the empirical value of 2.47 fm [44]. 

Figure 6 shows the charge form factors for the inter¬ 
band transitions between the two bands with n = 1 and 
n = 2 in Since the 0^ state is the Holye state with 
dilute 3 q: structure, which is beyond the model space of 
the present calculation, the inelastic form factor Fo{q) 
corresponding to the transition from the 0]'' to the 0^ 
states is significantly underestimated in the high-g region 
beyond the first minimum. It is worthwhile to mention 
that the calculated electric monopole transition matrix 
element \M{E0 : Oj —>■ 0]'')] = 4.1 efm^ and the charge 
radius of Oj (2.73 fm) are in good agreement with the 
results (4.5 ± 0.2 efm^ and 2.73 ± 0.02 fm, respectively) 
of the recent configuration mixing calculation based on 
a Skyrme force [43]. These values should be compared 
with the data \M{E0 : Oj —>■ 0]^)] = 5.4(2) efm^ and 
the charge radius from other calculations for the Holye 
state, such as 3.27 fm by the antisymmetrized molecular 
dynamics [45], 3.38 fm by the fermionic molecular dy¬ 
namics [46] and 3.83 fm by the alpha-condensation model 
[47]. 

The low-lying spectrum of with MR-CDFT cal¬ 
culation is shown in Figure 7, in which the results are 
compared with the experiment data [48, 49] as well as 
with other model calculations[43, 45, 46]. One can see 
that the low-lying spectrum is reproduced rather well by 
the present calculation, although the excitation energies 
are systematically overestimated. The electric monopole 
transition matrix element \M{E0)\ and the quadrupole 
transition strength B{E2) are also in good agreement 
with the data and the other model calculations. 
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r (fm) r (fm) 

FIG. 4: (Color online) The vector transition density, given by Eq. (13a), and the scalar transition density, given 

by Eq. (13b), in the low-lying states (n = 1) of These are plotted by the solid and the dashed lines, respectively. The 
insets show the difference of the vector and scalar transition densities. 


2. AN interaction 

The core states obtained in the previous subsection 
are used as inputs for the coupled-channels calculations 
for the hypernucleus To this end, the parameters 
and ay^ in the NA interaction, Eqs. (6) and 
(7), are fitted with the microscopic particle-rotor model 
to the experimental A binding energy of that is, 
^^exp.) _ 38 j-Q Q 5 MeV [1]. Figure 8(a) shows a con¬ 

tour plot of the absolute value of the difference between 
the theoretical and the experimental hyperon binding en¬ 
ergies as a function of and ay^. This is obtained 
by including in Eq. (1) the core states up to ricut = 2 
and /cut = 4. Obviously, the two strength parameters 
cannot be uniquely determined by fitting only to Ba and 
are linearly correlated as illustrated in Fig. 8 (a). 

Taking a few sets of the parameters along the valley 
with Bf = in Fig. 8(a), we calculate the energy 

of each of the low-lying excited states of 3/2+, 3/2“ and 
1/2“ in (see Fig. 8(b)). One can see that the exci¬ 
tation energies of 3/2+, 3/2“ depend on the choice of 
the parameters only weakly. The energy of 1/2“ state 
slightly decreases with the decrease of the absolute value 
of the coupling strengths. For all the sets of the pa¬ 
rameters {ag^,ay^) in the region of concerned, the en¬ 
ergy splitting between the first 1/2“ and 3/2“ states 
is in agreement with the data 152 ± 54(stat) ± 36(syst) 


keV [53, 54]. Therefore, as one of the choices, we first fix 
the value of to be —4.2377 x 10“® MeV“^, which is 
the same value as in the PCY-S2 set [28], and determine 
ay-^ to be 1.969 x 10“^ MeV“^. With this parameter set 
for the AN interaction, the energy splitting between the 
1/2“ and 3/2“ states is 198.9 keV. 


3. Projected potential energy surface 

With the NA interaction determined in the previous 
subsection, let us first investigate the projected potential 
energy surface for Figure 9 shows the resultant en¬ 
ergy Ej{j3) for the = 1/2+, 1/2“ and 3/2“ states in 
as a function of the deformation j3 of the core nucleus. 
These are obtained by solving the coupled-channels equa¬ 
tions for each /3. For comparison, the figure also shows 
the potential energy curve in the single-channel calcu¬ 
lations without taking into account the core excitations 
(the dot-dashed lines) as well as those for the core nu¬ 
cleus (the dashed lines). The energy surfaces obtained 
with the coupled-channel calculations are systematically 
lower than that with the single-channel calculation due 
to the additional configuration mixing effect. One can 
see that the hypernuclear energy curve with spin-parity 
of 1/2+ has an oblate minimum with j/3j significantly 
smaller than that of with 0+, indicating a smaller 






























q 

FIG. 5: (Color online) The charge form factor for the tran¬ 
sition from the ground state to (a) the ground state, (b) the 
excited 2f state, and (c) the excited 4^ state in calculated 
with the GCM method with the PC-Fl force, in comparison 
with the available data [50, 51]. R on the top panel is the 
root-mean-square charge radius of 



q (fm ’) 

FIG. 6: (Color online) The charge form factors between the 
two bands with n — I and n = 2 in ^^C calculated with the 
GCM method with the PC-Fl force, in comparison with the 
available data from Ref. [46]. 
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FIG. 7: (Golor online) The spectrum of ^^C obtained with 
several methods, that is, the GGM with the PC-Fl interac¬ 
tion (the present calculation), GCM with the Skyrme SLy4 
force [43], anti-symmetrized molecular dynamics (AMD) [45] 
and the alpha cluster model [46]. The excitation energies 
are given in units of MeV. The solid and the dashed arrows 
are the quadrupole transition strength B{E2) (e^ fm"*) and 
the electric monopole transition matrix element \M{E0)\ (e 
fm^), respectively. The experimental data are taken from 
Refs. [48, 49]. 



(10 " MeV") a"* (10 " MeV ") 


FIG. 8: (Color online) (a) A contour plot of the absolute 
value of the difference between the theoretical and the exper¬ 
imental hyperon binding energies of hypernucleus as a 
function of the coupling strength parameters and aY^ 
in the NA interaction. The calculated energies are obtained 
with the microscopic particle-rotor model, (b) Energy levels 
of the 3/2^, 3/2“ and 1/2“ states in C calculated with the 
strength parameters denoted by the dots in the panel (a). 


collectivity. On the other hand, for J'^ = 1/2“, the de¬ 
formation around the oblate minimum is similar to that 
of but with a higher barrier at the spherical shape. 
This leads to a smaller effect of shape mixing between 
the prolate and oblate configurations and thus a larger 
average deformation in The main component of the 
1/2+ and 1/2“ hypernuclear states are the A particle in 
Si /2 and pi /2 orbits coupled to the ground state of 
respectively. It implies that a A particle in the s (p) or¬ 
bit decreases (increases) the collectivity of which is 
consistent with the findings in the recent studies [14, 23]. 
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FIG. 9: (Color online)The potential energy surfaces of hyper¬ 
nucleus obtained with the single-channel calculation (the 
dot-dashed lines) and the coupled-channel calculation (the 
solid lines) with spin-parity of (a) = 1 /2’'', (b) = 1 /2“ 

and (c) = 3/2“ as a function of deformation parameter 

/3. For comparison, the figure also shows the energy surface 
for the nuclear core with spin-parity of = O’*' (the dashed 
lines). The energy surfaces for hypernuclei are shifted by a 
constant amount as indicated in each panel. 


4- Single-channel calculations 

Let us now employ the core states described with the 
GCM method and discuss the spectrum of Before we 
present the results of full coupled-channels calculations, 
we first discuss the results of single-channel calculations, 
restricting the A-hyperon to a specific orbital coupled to 
a single core state. In this case, we take only the diagonal 
component in the coupling potentials, Eqs. (11) and (12), 
in the coupled-channels equations. The coupled-channels 
equations are then simplified as 

9k{r) + {Ei^ - Ej)fk{r) 

+ [C/f(r) + C/f(r)]/,(r) = 0, (28) 

f fkir) - (E/„ - 2mA - Ej)gk{r) 

\dr ^ / 

-[C/f(r)-t/f(r)]gfe(r)=0. (29) 

The results for the A particle in the Si/ 2 , P 1/2 and 
P 3/2 orbitals are shown in the columns (d),(e),(f) and (g) 
of Figure 10. For comparison, the figure also shows the 
spectrum of the core nucleus in the columns (a), (b), 
and (c) (these are actually the same as those in Fig 7). A 
hyperon coupled to the core state with angular mo¬ 
mentum /"*" produces several hypernuclear states with 
with the total angular momentum J running from 
J = |/ — j| to J = / -I- j, with the parity of tt = (—1)^. 
When the A particle is restricted to the Si /2 orbit, a dou¬ 
blet states with (/ — 1/2)+ and (/ -|- 1/2)+ are yielded, 
which are degenerate in energy for / > 0. On the other 


hand, a spectrum is more complex for the case of A par¬ 
ticle in the P 3/2 orbital. In this case, the multiplet states 
with J € [\I — j\, {I 4- j)] are ordered according to the 
properties of the coupling potential in Fq. (28). 

For instance, in the case of Ap 3/2 ® 2+, the multiplets 
are ordered as 5/2“, 3/2“, 7/2“ and 1/2“ (see the col¬ 
umn (g) in Fig. 10). In order to understand this, we 
write the coupling potentials for the configuration k as 

[yfcfc = (-i)^+^+"^(jT||y.||jT)|;( ^ 

= (30) 

A 

where the indices m = S and V represent the scalar and 
vector potentials, respectively. The coefficients Cn {n = 
0,1, and 2) are defined as Cq = (—1)1+^+'^, 


Cl A = (j^||>A||j^) 
(-1)^ + V2 


\/47r 




(j A 

1 ^ 1/2 0 


j 

- 1/2 


and 


C2A = 


J I J 
A J I 



(31) 

(32) 


Table I lists the value of each of the coefficients Cn- 
The transition densities p^^^(r) and the potential 
U^^{r) + Uf^{r) with k = = (|,l,2i) are dis¬ 

played in Fig. 11 as a function of radial coordinate r. It is 
seen that the potential Uy^{r) -\- Ug^{r) becomes gradu¬ 
ally deeper in the order of J'^ = l/2“, 7/2“, 3/2“, 5/2“, 
which is consistent with the distribution of energy lev¬ 
els of these multiplets. It is seen in Table I that the 
product C 0 C 1 AC 2 A is the same among the multiplets for 
A = 0, and thus the origin for the energy difference among 
these four hypernuclear states is the non-zero A = 2 
term in the potential (30). That is, the splitting of 
= l/2“, 7/2“, 3/2“, and 5/2“ hypernuclear states 
is originated from the non-zero transition density p^{r) 
(cf. Fig. 11(a)) due to the reorientation effect (that is, 
the transition between the same state) of 2+ state in the 
deformed shape of Since negative as 

shown in Fig. 11(a), the potential is most attractive for 
J'^ = 5/2“, which has a positive C 0 C 12 C 22 . 

For the configuration with A in Pi /2 orbital coupled 
to the nuclear core 2+ state, the resultant doublet states 
3/2“ and 5/2“ are degenerate in energy, since the coef¬ 
ficient C = C 0 C 1 AC 2 A is not zero only for A = 0, having 
the same value of 1 /a/ 20^ between those two states (see 
Table I). 


5. Coupled-channels calculations 

Let us now solve the coupled-channels equations for 
^aC. To this end, we first examine the convergence fea¬ 
ture of the excitation energies with respect to the cutoff 
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FIG. 10: (Color online) The low-energy excitation spectra of ((a)-(c)) and ((d)-(l)). For the full GCM calculations 
are compared with the experimental data taken from Ref. [52]. The columns (h) and (i) show the positive-parity states in ^aC, 
while the columns (j) and (k) show the negative-parity states. The experimental data of ^aC are taken from Refs. [53, 54]. 
For comparison, the results of single-channel calculation for ^aC with the A particle in the Si/ 2 , P 1/2 and P 3/2 orbitals are also 
plotted in the columns (d), (e), (f), and (g). 


TABLE I: The coefficients in the potential for the [Ajj 02"'"]^'^^ 
configurations (see Eqs. (31) and (32)). 


[Aij (g) 2+] r 

Go G 20 G 0 G 10 G 20 

G 22 G 0 G 12 G 22 

[Ap3/2(g)2+] 1/2“ 

[Ap3/2 <8) 2+] 3/2“ 
[Ap3/2 » 2+] 5/2“ 
[Ap3/2 ® 2+] 7/2“ 

1-00 v4o 

1-00 v4o 

1-00 io 

1-00 ,Ao 

20 20 

0.00 0.00 

vT4 

28 28 

70 70 

[Api/2 g) 2+] 3/2“ 

[Api/2 ® 2+] 5/2“ 

1-00 io 

,/Tn 

0.00 0.00 

0.00 0.00 


of core states ricut and the core angular momentum /cut- 
Fig. 12 shows that ricut = 2 and /cut = 4 yield a good 
convergence for the low-lying excited states, and we use 
these cut-offs in the calculations presented below. 

The columns (h), (i), (j), and (k) in Figure 10 show 
the calculated low-energy excitation spectra of in 
comparison with the corresponding data. One can see 
that the low-lying spectra for are reproduced rather 
well, although the excitation energies are slightly overes¬ 
timated. 

In the coupled-channels calculation, the doublets 
(5/2+, 3/2+) and (9/2+,7/2+) in the column (h) mainly 
consist of the configuration of Asi/ 2 ® 2 + and Asi/2(8)4+, 
respectively. See Tab. II for the probabilities for the 
dominant components in each state. These doublets are 



r (fm) r (fm) 

FIG. 11: (a) The transition densities ovPaV "f (-^ = 0 

and A = 2) and (b) the potential Uy*‘{r) + Ug'‘{r) in Eq. 
(28) for the hypernuclear states J = 1/2“, 3/2“, 5/2“, 7/2“ 
with the Ap 3/2 O 2+ configuration as a function of the radial 
coordinate r . 


degenerate in the single-channel calculation, as already 
shown in the column (d) in Fig. 10. The states of 3/2+ 
and 5/2+ are different from each other by 10 keV due to 
the weak mixing of other configurations. The levels in 
the column (i) correspond to the configuration of Asj ^/2 
coupled to the second band (n = 2) in These levels 
share similar features as those in the column (h). 

The negative-parity states are shown in the columns (j) 
and (k) in Fig. 10. The energy splitting between the 3/2“ 
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cut cut 

FIG. 12: (Color online) Excitation energy of low-lying states 
in as a function of the cutoff of core states n ((a) and 
(b)) and the cutoff of core angular momentum 7 ((c) and (d)) 
for the coupled-channels calculations. 


and 1/2“ states is as small as 199 keV. Notice that in the 
single-channel calculation the energy difference between 
the pure configurations of Ap 3/2 G) O/" and Api /2 ® 0 )^ 
is 180 keV. That is, the energy splitting of 3/2“ and 
1/2“ states reflects mainly the spin-orbit splitting of A 
hyperon in the P 3/2 and pi /2 states. A small splitting 
between the 3/2“ and 1/2“ states has been shown also 
in our previous calculation for ® Be [27], although it does 
not reflect the Lambda spin-orbit splitting because of a 
strong mixing between the Ap^^^ ® 0 ^ and the Apg^^ ® 2 + 
configurations in the 1 / 2 “ state. 

For the second 1/2“ and 3/2“ states, one can see a 
large configuration mixing (see Table II). This is because 
there are two states whose unperturbed energy in the 
single-channel calculations, close to one another. 

These two states are strongly coupled due to the off- 
diagonal components of the coupling potentials in the 
coupled-channels equations. Notice that, in this 
happens already in the first 1/2“ state [27], because the 
reorientation effect discussed in the previous subsection 
brings the Apg^^ ® 2 + configuration close to the Ap^^^ ^ 0 '*' 
configuration in energy due to the prolate nature of the 
2+ state of ®Be. 

According to our calculation, the experimentally ob¬ 
served level at excitation energy of 11.8 MeV has the 
spin-parity of 1 / 2 “'", dominated by the configuration 
Asi /2 ® 0^ (c.f. Fig. 10(e)) or the first radial excitation 
state of the configuration kg-Yji ® 0/“ (c.f. Fig. 10(d)). 

Figure 13 shows a comparison of low-energy excita¬ 
tion spectra of obtained with the present micro¬ 
scopic particle-rotor model (MPRM), the multi-channel 
algebraic scattering (MCAS) approach [55], and the 
3a -b A cluster model [56] , together with the experimen- 


TABLE II: The probability Pjii^ for the dominant compo¬ 
nents in the wave function for low-lying states of ^iC ob¬ 
tained by the microscopic particle-rotor model. Only those 
components which have Pjlln larger than 0.1 are shown. E 
is the energy of each state obtained by solving the coupled- 
channels equations, while is the unperturbed energy ob¬ 
tained with the single-channel calculations. The energies are 
listed in units of MeV. 
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tal data [1]. The basic idea of MCAS approach for A 
hypernuclei[55] is similar to the microscopic PRM model 
in a sense that the hypernuclear wave function is given 
by the A hyperon coupled to the low-lying states of nu¬ 
clear core. In contrast to our full microscopic models, in 
which all the inputs are from the multi-reference CDFT 
calculation, the MCAS approach adopts the experimen¬ 
tal data for the energies of nuclear core states with an 
assumption of a pure collective rotational states and a 
phenomenological deformed Woods-Saxon potential for 
the coupling potentials. In contrast to the MCAS ap¬ 
proach and cluster model calculation, the ordering of the 
first degenerate 3/2“*“ and 5/2+ states are not reproduced 
in the microscopic PRM calculation, since we do not in¬ 
clude a spin-spin interaction in our calculation. Except 
for this, the ordering of low-lying states and the structure 
of spectrum are the same between the microscopic PRM 
and the MCAS approach. The main components of each 
state obtained with the microscopic PRM calcidation are 
similar to those in the cluster model calculation. 

Table III shows the calculated E2 transition strengths 
for low-lying positive parity states of the hypernucleus 
and the corresponding core nucleus. In order to remove 
the trivial factor due to the angular momentum coupling 
for Si /2 for the A particle and see more clearly the impu¬ 
rity effect of A particle on nuclear collectivity, we define 
the cB{E2) value (that is, the B{E2) value for the core 
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FIG. 13: (Color online) A comparison of low-energy excitation 
spectra of obtained with the present microscopic particle- 
rotor model (MPRM), the multi-channel algebraic scattering 
(MCAS) approach [55], the 3 q-|-A cluster model [56], and the 
experimental data [1]. 


TABLE III: The calculated E2 transition strengths (in units 
of fm'*) for low-lying positive parity states of and ^iC. 
The cB{E2) values are calculated according to Eq. (33). The 
changes in the B{E2) is indicated with the quantity defined 
by A = {cB{E2) - B{E2-^^C))/B{E2-,^^C). The value in 
the parenthesis for is the experimental data taken from 
Ref. [48]. 
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9/2+ 
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11.48 

11.48 

-21.36 


part) as, 


cB{E2 : li If) 

= ^2*1 B{E2 : J, ^ Jf), (33) 


where ji is the value for the main channel in the initial 
state. The impurity effect of A particle can be discussed 
by comparing the B{E2) values for the core nucleus and 
the cB{E2) values for the corresponding hypernucleus. 
One can see that the E2 transition strength for 2^ —0]'' 
in is significantly reduced, by a factor of ^ 14%, due 
to the addition of a A particle. 



Exp. PC-F1 PC-PK1 Cluster AMD 


FIG. 14: (Color online) A comparison of the yrast rotational 
states of ^°Ne obtained with several methods. The results of 
the cluster model and the AMD are taken from Refs. [57] 
and [14], respectively. 


B. Low-energy spectroscopy of "iNe 

We next consider an application to hypernuclei in the 
sd-shell region. For this purpose, we discuss the ^^Ne hy¬ 
pernucleus. Since the A binding energy in has not 
yet been measured, we fit the value of coupling strength 
parameters and to the A binding energy es¬ 
timated with a deformed relativistic mean filed calcula¬ 
tion [23]. With the PC-Fl and PCY-Sl forces for NN 
and NA interactions, respectively, i?A is estimated to be 
14.35 MeV for the lowest A hyperon state. With the same 
process as in Sec. Ill A-2, we obtain a parameter set of 
= -4.2377 X 10-5 MeV^ and = 1.6694 x lO"® 
MeV-2. 

Figure 14 shows the calculated yrast rotational states 
of ^°Ne. In order to see the parameter set dependence 
for the NN interaction, we use both PC-Fl and PC-PKl 
parameter sets. For a comparison, the figure also shows 
the results of a+^^0 cluster model [57] and the AMD 
model [14]. One sees that all these models reproduce the 
rotational character of the yrast states, although they 
tend to overestimate the moment of inertia. 

Figure 15 shows the obtained energy curve Ej{f3) for 
the = 1/2“*“ and 1/2“ states in ^^Ne as a function of 
the deformation /3 of the core nucleus. The left and the 
right panels show the result with PC-Fl and PC-PKl 
forces, respectively. For the latter, we use the same NA 
interaction as in the former calculation, even though the 
parameters are determined with PC-Fl. We have con¬ 
firmed that this yields the B\ value of 14.33 MeV with 
PC-PKl, which is similar to the value with the PC-Fl 
set, that is, 14.35 MeV. For PC-Fl, the hypernuclear 
energy curve with spin-parity of 1/2+ and 1/2“ has a 
prolate minimum with a smaller /3 than that of ^°Ne 
with 0+. For PC-PKl, on the other hand, the value 
of j3 at the energy minimum remains almost the same 
for the 1/2+ conhguration while that for the 1/2“ con- 
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FIG. 15: (Color online) Energy curve Ej{/3) for the J'^ = 
1/2+ (the solid line) and 1/2“ (the dot-dashed line) states in 
^^Ne as a function of the deformation /3 of the core nucleus. 
These are obtained with PC-Fl (the left panel) and PC-PKl 
(the right panel) forces. In order to make a comparison easy, 
each hypernuclear curve is shifted by a constant value so that 
the energy at the absolute minimum coincides with that for 
^°Ne with 7'^ =0+. The difference between the energy curve 
of ^aNo and that of ^°Ne is shown in the panels (c) and (d) 
for the PC-Fl and PC-PKl forces, respectively. 


figuration increases as compared to the deformation for 
^°Ne with 0+. Notice that the energy surface for 1/2“ 
has a higher barrier at the spherical shape than the bar¬ 
rier for ^°Ne for both the interactions. This indicates 
that with 1/2+ and 1/2“ has a smaller and a larger 
collectivity than that of ^°Ne, as in The energy dif¬ 
ferences between the 1/2+ state in ^^e ground 

state of ^°Ne, as well as the 1/2“ state in ^j^Ne and the 
ground state of ^°Ne, are shown in Figs. 15 (c) and 
15 (d). Even though PC-Fl and PC-PKl forces predict 
somewhat different energy curves, those energy curves 
are qualitatively similar to each other, especially when 
those are plotted with respect to the energy curve for 
20Ne. 

Figure 16 shows the spectra of both ^°Ne and ^iNe. 
To examine the channel-coupling effect on hypernuclear 
states, we also include the results from single-channel cal¬ 
culations. The figure only shows the results with PC-Fl, 
since the results with PC-PKl are similar (see also Fig. 
18 below). The probability of the main components is 
summarized in Table IV. It is shown that the hypernu¬ 
clear states with positive parity plotted in the columns 
(i) and (j) from the full coupled-channels calculation are 
close to the results of single-channel calculation shown in 


TABLE IV: Same as Table II, but for ^^Ne hypernucleus with 
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the columns (e) and (f). The analysis of hypernuclear 
wave functions demonstrates that these states are domi¬ 
nated by the configuration of /2 coupled to the state I 
in the first {n = 1) and the second (n = 2) bands in ^°Ne, 
respectively, with the weight between 98% and 99%. It is 
seen that the hypernuclear doublet states 
with configuration A^i/ 2 C)d+ are degenerate, as discussed 
in Sec. Ill A 4. Moreover, it is seen that the spectra of 
positive-parity states in ^a^c is close to that of ^°Ne with 
similar excitation energies to each other. In other words, 
the presence of a Asi /2 does not change significantly the 
low-energy structure of the core nucleus ^°Ne. 

The negative-parity states in ^ANe are shown in the 
columns (k) and (1). One can see that the channel¬ 
coupling effect plays an important role in their excitation 
energies. Moreover, we note that the energy difference 
between the first 1/2“ and 3/2“ states is less than 40 
keV. Notice that the 1/2“ state is a strong admixture of 
the configurations Api /2 ® 0)*' and Ap 3/2 0 On the 
other hand, the 3/2“ state is a strong admixture of the 
configurations Ap 3/2 0 0)*", Ap 3/2 0 2+ and Api /2 ® 2+. 
Therefore, the splitting of the 1/2“ and 3/2“ levels in 
^ANe does not reflect the strength of A spin-orbit inter¬ 
action, which is in marked difference from the case in 
^aC. From yet another point of view, it is interesting to 
point out that a typical rotational band having L=l“, 
3“, 5“, • • • is realized as seen in the column (k) of Fig. 
16, apart from the spin of the hyperon. This group can 
be charactrized by the K = 0~ band based on the strong 
coupling between the nuclear rotation and the hyperon 
in the p-state, and thus this band manifests a genuinely 
hypernuclear state with the [5] (90) symmetry which is 
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FIG. 16: (Color online) The low-energy excitation spectra of ^°Ne and ^iNe obtained with the microscopic particle-rotor model 
calculations. 



FIG. 17: (Color online)Same as Figure 11, but for ^i^e. 


similar to the [5] (50) band verified in ®Be [3]. This fea¬ 
ture of negative parity states in hypernuclei has already 
been discussed in Ref. [58] with the traditional particle- 
rotor model with the Elliot SU(3) model for the core 
states. 

The underlying reason for the difference between 
and ^o the different properties of the core 

nuclei. ^®Ne is well-deformed with a much larger tran¬ 
sition density P 2 ^{r) than that in Notice also that 
the ordering of the Ap 3/2 0 2'^ multiplet states is op¬ 
posite to that in reflecting the fact that the sign 
of quadrupole moment is opposite (that is, prolate de¬ 
formation for ^°Ne and oblate deformation for ^^C). In 
Figure 17, we plot the transition density for the ^iNe 
hypernucleus. One can see that the A=2 component has 
the opposite sign as compared to the transition density 


for shown in Fig. 11. Similar to the ®Be case, these 
result in several 1/2“ and 3/2“ states close in energy in 
the single-channel calculations, which are strongly mixed 
in the full coupled-channels treatment. Similarly to the 
1/2“ and 3/2“ states, the 5/2“ and 7/2“ states also 
show the strong configuration mixing between A-hyperon 
in Pi /2 and P 3/2 orbits. We have found that this feature 
of strong mixing found in ® Be and persists also in 
heavier systems, such as 31 Si and is^Sm. 

Notice that the ^°Ne nucleus has prominent negative- 
parity bands originated from the q;-|-1®0 structure. For 
simplicity, in the present calculations, we have assumed 
reflection symmetry for ^°Ne. The inclusion of these neg¬ 
ative parity states in the coupled-channels calculations is 
thus beyond the scope of the present paper. It would 
be an interesting future work to include them and study 
how the negative parity states in Are perturbed. 

Figure 18 shows a comparison of low-energy excitation 
spectra of ^^Ne obtained with the cluster model [57], the 
AMD [14], and the present microscopic PRM calcula¬ 
tions based on the PC-Fl and PC-PKl interactions. The 
positive-parity band in the microscopic PRM is closer to 
the result of cluster model as compared to the result of 
AMD, which has a slightly larger moment of inertia. The 
negative-parity states are similar to the AMD results but 
with lower excitation energies, which might be due to the 
large channel-coupling effect taken explicitly into account 
in the present work, see Fig. 16 and Table IV. 

Table V lists the E2 transition strengths for low-lying 
states of with PC-Fl. For comparison, the table 

also shows the change in B{E2) from ^°Ne to ob¬ 
tained with the PC-PKl force, the cluster model [57], and 
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FIG. 18: (Color online) A comparison of low-energy excita¬ 
tion spectra of obtained with the cluster model [57], the 
AMD [14], and the present microscopic particle-rotor model 
(MPRM) calculations with the PC-Fl and PC-PKl forces. 


TABLE V: The calculated E2 transition strengths (in units of 
fm^) for low-lying states of ^l^e with the PC-Fl force for 
the core states. The results for the change in the B{E2) value 
from ^°Ne to ^iNe is compared with the results with PC-PKl, 
the cluster model [57] and the AMD [14] calculations, where 
A is defined in the caption of Table III. 


Transition 


PC-Fl 


PC-PKl AMD Cluster 


3/^ 

5/2+ 

7/2+ 

7/2+ 

9/2+ 


Jf JJ B{E2) cB{E2) A(%) A(%) A(%) A(%) 


1/27 

1 / 2 + 

3/2+ 

5/2+ 

5/2+ 


54.28 

54.28 

-3.19 

-7.16 

-11.8 

-23.9 

54.28 

54.28 

-3.19 

-7.16 

-11.5 


65.90 

73.22 

-3.95 

-4.80 

-17.8 

-22.6 

7.32 

73.22 

-3.95 

-4.80 



73.22 

73.22 

-3.95 

-4.81 

-13.0 



the AMD [14]. The B{E2) value decreases by adding 
a A hyperon in s-orbit in these calculations. However, 
the cluster model and AMD model predict more reduc¬ 
tion compared to the microscopic PRM. A further study 
would be necessary in order to reconcile this difference. 

We have also applied the microscopic particle-rotor 
model to study another sd- shell hypernucleus, ^Si. We 
have found that the impurity effect of A hyperon in 7si 
is qualitatively the same as that in 7Ne. 


C. Low-energy spectroscopy of Sm 

One of the advantages of the microscopic particle-rotor 
model is that this method is not limited to light hy¬ 
pernuclei but it can also be applied to medium-heavy 
and heavy hypernuclei. As an example of application to 


heavy deformed hypernuclei, we next consider ^®^Sm and 
^^^Sm. By fitting to Ha=24.98 MeV estimated with the 
deformed RMF calculation, we obtain a parameter set of 
AN interaction as ag^ = —4.2377 x 10“® MeV“^ and 
ay^ = 1.0401 X 10“®MeV~^. Figure 19 shows the pro¬ 
jected energy curves for ^®^Sm and ^^^Sm obtained with 
this NA interaction together with PC-Fl for the NN in¬ 
teraction. For the 1/2+ state, the polarization effect of 
A particle in s-orbit on the properties of ^^^Sm is much 
smaller than that on and ^°Ne due to the large mass 
number, although the effect is still large for the negative 
parity states due to the strong channel coupling effects. 

Figure 20 shows the calculated low-energy spectrum of 
^®^Sm and ^®^Sm with the PC-Fl force. The ground- 
state band and the two /3-bands in ^^+Sm are rea¬ 
sonably reproduced, although the band-head energy of 
the /3-bands are overestimated. The low-lying positive- 
parity states in ^^^Sm are dominated by the single¬ 
configuration of Asj ^2 0/+ with similar excitation energy 
as that of the nuclear core state with /+. As shown in 
Tab. VI, the positive-parity states J+, except for 1/2+, 
are nearly two-fold degenerate. These characters are sim¬ 
ilar to the hypernuclei in the light-mass region. On the 
other hand, one can see that the negative-parity bands 
are well separated from the positive-parity ground band 
in ^®^Sm, which is different from the light hypernuclei. 
It is because the energy scale of the rotational motion 
is proportional to (see Eq. (1.50) in Ref. [38]), 

while that of single-A excitation from s to p orbit is pro¬ 
portional to . Therefore, with the increase of mass 
number A, the rotational energy spectrum is compressed 
faster than the single-A excitation spectrum. Besides, the 
low-lying negative-parity states J~ are nearly two-fold 
degenerate, even though there are strong configuration¬ 
mixing in these states. 

Table VII presents the E2 transition strengths in 
^®^Sm and ^®^Sm. It is shown that the change in the 
B{E2) values by adding a A hyperon in s orbital is less 
than 1%, which is much smaller than that in the light hy¬ 
pernuclei studied in this paper. This is consistent with 
the small polarization effect of A particle discussed in 
connection to the projected energy surface shown in Fig. 
19. 


IV. SUMMARY 

We have presented the detailed formalism of the mi¬ 
croscopic particle rotor model based on a covariant den¬ 
sity functional theory for the low-lying states of single-A 
hypernuclei. In this formalism, the wave functions for 
hypernuclei have been constructed by coupling the A hy¬ 
peron to the low-lying states of the core nucleus. The 
radial wave functions are obtained by solving the corre¬ 
sponding coupled-channel equations, in which the cou¬ 
pling potentials are provided in terms of the transition 
densities of the nuclear core states. For simplicity, in 
this paper we have adopted only the leading-order four- 
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FIG. 19: (Color onlme)Same as Figure 9, but for and 



FIG. 20: (Color online) The low-energy excitation spectra of 
((a)-(f)) and «lSm ((g)-(i)). 


fermion coupling terms of scalar and vector types for the 
AN effective interaction. Applying this method to 
we have reproduced reasonably well the experimental en¬ 
ergy spectrum of this hypernucleus. We have applied this 
method also to ^^^e, and and have achieved a 

good agreement with other model studies both for the ex¬ 
citation energies and the compositions of wave functions. 
We mention that our method is the only one which can 
be applied to such heavy hypernuclei. We have found 


TABLE VI: Same as Table II, but for ^^ISm hypernuclei. 


r 

E 

(O') 


r 

Pjlln 

^( 0 ) 

-^Ich 

1/2+ 

0.00 

Sl/2 


ot 

0.98 

0.00 

^72^ 

0.11 

Sl/2 



0.98 

0.11 

5 / 2 + 

0.11 

Sl/2 


2/ 

0.98 

0.11 


0.35 

Sl/2 


4 

0.98 

0.35 

9 / 2 + 

0.35 

■81/2 


4 

0.98 

0.35 


1.68 

81/2 


02+ 

0.99 

1.58 


1.99 

81/2 

® 

4 

0.99 

1.90 

5/22+ 

1.99 

81/2 

® 

4 

0.99 

1.90 

1/27 

6.28 

P3/2 

® 

4 

0.66 

7.58 



Pl/2 

® 

0/ 

0.32 

8.51 

3/27 

6.27 

P3/2 

® 

0/ 

0.35 

8.31 



P3/2 

® 

2/ 

0.32 

8.42 



Pl/2 

® 

2/ 

0.29 

8.61 

5/27 

6.44 

P3/2 

® 

4 

0.53 

8.08 



Pl/2 

® 

4 

0.33 

8.61 



P3/2 

® 

4 

0.11 

8.96 

7/27 

6.44 

P3/2 

® 

4 

0.47 

8.19 



Pl/2 

® 

4 

0.28 

8.86 



P3/2 

® 

4 

0.22 

8.89 


TABLE VII: Same as Table III, but for ^^^^Sm and i®|Sm. 
The experimental data for ^®^Sm, shown in the parenthesis, 
is taken from Ref. [59]. 


^^^Sm 

If If B{E2) 

rsm 

Jf Jf B{E2) cB{E2) A(%) 

2+ 0/ 9358.49 

3/2+ 1/2+ 9284.69 9284.69 -0.79 

(8720 ± 100) 

5/2+ 1/2+ 9284.23 9284.23 -0.79 

4/ ->■ 2/ 13512.18 

7/2/ 3/2+ 12081.04 13423.38 -0.66 


7/2/ 5/2/ 1342.25 13422.54 -0.66 


9/2+ ->■ 5/2+ 13422.25 13422.25 -0.67 


that the NN interaction with the PC-Fl and PC-PKl 
sets lead to similar hypernuclei spectra to each other. 
For all the hypernuclei, the low-lying excited states with 
positive parity J'*', except for 1/2’'', are nearly two-fold 
degenerate and dominated by the single-configuration of 
As ,^^2 ® where /“*■ is the spin-parity of the nuclear 
core states. In contrast, in general there are large con¬ 
figuration mixing in the negative-parity states. We have, 
however, found an exception for this, that is, for the first 
3/2“ and 1/2“ states in the effect of configuration 
mixing is rather small, and thus the energy splitting of 
these states reflects the spin-orbit splitting of A hyperon 
in the P3/2 and pi/2 states. Concerning the electromag¬ 
netic transitions, we have found that for all the systems 
the B{E2) value from the first 2+ to the ground states in 
the core nuclei is reduced by adding a A particle in the 
positive-parity states. The reduction factor is about 14% 
for 3.2% for ^/^Ne, and 0.79% for ^®^Sm, and thus 
the reduction factor is larger for the oblate hypernuclei. 
For ^ slightly larger impurity effect was 
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found with the PC-PKl force as compared to the PC-Fl 
force. 

New measurements of 7-ray spectroscopy of hypernu¬ 
clei will soon start at the new generation experimental 
facilities such as J-PARC. It would be interesting if the 
low-lying spectra predicted in this paper are confirmed 
in near future. 
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